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Abstract 

We apply perturbative QCD factorization theorems to inclusive heavy 
hadron decays, and obtain simultaneously a low semileptonic branching ratio 
B(B — > Xlv) = 10.16%, the average charm yield (n c ) = 1.17 per B decay, 
a small lifetime ratio r(Ab)/r(B c i) = 0.78, and the correct absolute decay 
widths of the B meson and of the At, baryon. 
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Dynamics of heavy hadrons is greatly simplified in the limit of infinite 
quark mass. Heavy quark effective theory (HQET) provides a systematic ex- 
pansion of QCD in this limit, which makes possible an almost-first-principle 
analysis of inclusive heavy hadron decays. Two different QCD-based ap- 
proaches have been developed to study these decays: the HQET-based op- 
erator product expansion (OPE) J]J and perturbative QCD (PQCD) factor- 
ization theorems |J. In the OPE formalism the relevant hadronic matrix 
elements are expanded in inverse powers of the heavy quark mass Mq. By 
contrast, in the PQCD formalism the factorization of a decay rate into a 
hard subamplitude and a universal heavy hadron distribution function is 
performed according to inverse powers of the heavy hadron mass Mh- It has 
been shown that to order 1/Mq the two approaches give identical results in 
the case of inclusive semileptonic decays [§]. 

One of the puzzles in heavy hadron decays is the explanation of the low 
lifetime ratio r(A 6 )/r (B d ) = 0.79±0.06 @. The OPE prediction to 0(1/M fe 2 ), 
Mfc being the b quark mass, is about 0.99 Q. When including the 0(1/Mjf) 
corrections, the ratio depends on six unknown parameters, and reduces to 
around 0.97 for various model estimates, which is still far beyond the cor- 
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rect value. It is known that the replacement of the overall M% factor in 
front of nonleptonic decay widths by the corresponding hadron mass Mfj 
provides a possible solution to the puzzle ||, since (M B /M\ b ) 5 = 0.73, M B 
and M^ b being the B meson mass and the A& baryon mass, respectively, is 
close to the observed ratio. As pointed out in ||, the same replacement also 
explains the absolute B meson decay rate, while the OPE approach using 
the expansion parameter M& accounts for only 80% of the decay rate. How- 
ever, the kinematic replacement for nonleptonic decays is at most an ad hoc 
phenomenological ansatz without solid theoretical base. 

In Ref. we have performed the PQCD factorization of inclusive heavy 
hadron decays according to hadronic kinematics, such that the above re- 
placement is implemented naturally. Generally speaking, a decay rate is 
expressed as the convolution of a hard b quark decay subamplitude with a 
universal heavy hadron distribution function and several jet functions for 
light energetic final states. A jet function contains double logarithms from 
the overlap of collinear and soft enhancements in radiative corrections. For 
the semileptonic decays b — > cW, because the final-state c quark is massive, 
there are no double logarithms and thus no jet functions. The transverse 
momentum p± carried by the b quark is introduced as a factorization scale, 
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above which PQCD is reliable, and radiative corrections are absorbed into 
the hard subamplitude, and below which physics is regarded as being purely 
nonperturbative, and absorbed into the distribution function. To facilitate 
the factorization, we work in the impact parameter b space, which is the 
Fourier conjugate variable of p±. 

For the nonleptonic modes b — > cud and b — > ccs, which involve more 
complicated strong interactions, additional subprocesses need to be included 
except for those in the semileptonic cases. These modes involve three scales: 
the W boson mass Mw, at which the matching condition of the effective 
Hamiltonian for nonleptonic decays to the full Hamiltonian is defined, the 
characteristic scale t of the hard subamplitude, which reflects the specific dy- 
namics of heavy hadron decays, and the b quark transverse momentum p±, 
or equivalently 1/6, which serves as the factorization scale mentioned above. 
To analyze the nonleptonic decays, we have proposed a three-scale factoriza- 
tion formalism that embodies both effective field theory and the standard 
PQCD factorization theorems. The hard gluon exchanges among quarks, 
which generate logarithms of Mw, are factorized into a harder function [^]. 
Collinear and soft divergences exist in radiative corrections simultaneously, 
when the light final-state quarks, such as u, d, and s, become energetic. The 
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resultant double logarithms ln 2 (p6), p being the jet momentum defined later, 
are absorbed into the corresponding jet functions. 

The various logarithmic corrections are summed to all orders using the 
resummation technique and renormalization-group (RG) equations, whose 
results are the evolutions among the scales Mw, t, and 1/6. The summation 
of the logarithms \n(Mw/t) is identified as the Wilson coefficient of the effec- 
tive Hamiltonian, which describes the evolution from the characteristic scale 
Mw of the harder function to the characteristic scale t of the hard subam- 
plitude. The summation of the logarithms \n(tb) leads to the evolution from 
t to 1/6. The double logarithms \n 2 (pb) in the jet functions are resummed to 
give a Sudakov factor 0, which suppresses the long-distance contributions 
from the large b region, and thus improves the applicability of PQCD. 

With the above ingredients, we have calculated the branching ratios of 
the modes b — > cW, b — > cud and b — > ccs for the inclusive B meson decays ||. 
It was observed that the decay rates of the b — ► cud and b —>■ ccs modes are 
both enhanced by the QCD evolution effects. In this way the semileptonic 
branching ratio B$l = B{B — > Xlv) is suppressed without increasing the 
average charm yield (n c ) per B decay. The controversy of the small B^l in 
inclusive B meson decays, which cannot be explained by the quark model, is 
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then resolved. This solution differs from the usual attempts in the literature, 

in which the b — > ccs mode is increased. In this letter we shall further show 

that our formalism provides a possible explanation of the low lifetime ratio 

r(A b )/r(B c i). Our predictions for the absolute decay widths of the B meson 

and of the baryon are also consistent with the data. 

We quote the factorization formulas for the semileptonic and nonleptonic 

B meson decays ||. The various momenta involved in the semileptonic 

decays B(P B ) — > X c + l(pi) + u{p u ) are expressed, in terms of light-cone 

coordinates, as 

( M B M B \ . , . , 

Pb = I -j= , -j= , Oj. I , Pi = {pi,Pi,0±), Vu = {Pv,P v ,Pv±), (1) 

where the minus component p~[ vanishes for massless leptons. For conve- 
nience, we adopt the scaling variables, 

_ 2E t _ q 2 _ 2q 

X ~W B " y -M%' m ~W B > {2} 

with the kinematic ranges, 

2y/a < x <l + a-P, 



,. _ N x + Vx 2 — 4a 

d < y < a + {1 + a — p — x) 



2 — x — y/x 2 — 4a 

x+ 2{V nA r^ y° <i+3/-A (3) 

x + v x z — Aa 
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where E% is the lepton energy and q = p\ + p u the momentum of the lepton 
pair. The constants a and (3 are a = Mf/M% and /3 = M^/M%, and M D 
being the lepton mass and the D meson mass, respectively. appears as 
the minimal invariant mass of the decay product X c . 

The b quark momentum is parametrized as Pb = Pb — P, where p = 
((1 — z)Pg , 0, p_i_) is the momentum of the light cloud in the B meson. The 
factorization formula for the total inclusive semileptonic decay width in the 
b space is given by || 

— = — - / dxdydy / dz dbbf B (z)J c (x,y,y ,z,b)H(x,y,y ,z) , 

1 AlX J JZmhx JO 

(4) 

with r = (G 2 F /lQiT 3 )\V c b\ 2 Ml. The momentum fraction z approaches 1 as 
the b quark carries all the B meson momentum in the plus direction. The 
minimum of z, determined by the condition p 2 c > M 2 , p c and M c being the c 
quark momentum and mass, respectively, is 



yo , m 2 x 



-y+ ' 



2 * M 2 Vx 2 -4a 

^min 



_yo + + ® 

2 x x 



1 yo x 



yo + y + ® 

2 x x 



(5) 



2 ^x 2 - 4a 

The function J c denotes the final-state cut on the c quark line, which is in 
fact part of the hard subamplitude. J c and the lowet-order H in momentum 



space are expressed as 



J c = S M B L-(l + z)^ + y + X{1 ~ Z) 



2 ^/x 2 - 4a 



2 x x 



-Ml - pi) , 



(6) 



X 



yx^ — Aa 



x \2\ l + z + {l - z) x — )-y- a )- ( ? ) 

The universal B meson distribution function fs, determined from the B — > 
Xs'j decay, is given by @, 

0.02647^(1 - zf 
JBy ) [(^- 0.95) 2 + 0.0034^] 2 ' K ) 

which minimizes the model dependence of our predictions. We have ignored 
the intrinsic b dependence of Jb from the infrared renormalon contribution, 
since it does not affect our predictions of the decay branching ratios. Note 
that the distribution function collects all-order soft gluons, which, being long- 
distance, are insensitive to the kinematics of the quarks they attach, such as 
masses. Hence, fs for the b — > 57 and for b — > c decays are the same. 
Through the z dependence of fs, the soft spectator effect is taken into ac- 
count in our formalism. 



Ignoring the penguin operators, the effective Hamiltonian for the nonlep- 
tonic decay b — > ccs is 

AC 

H cS = W c * s [ ci (m) Oi (m) + c 2 ( M ) 2 ( M ) ] , (9) 

with Oi = (s L r r^b L )(c L j' J 'CL) and 2 = {cLl f ib L ){sLl' 1 c L ). For the 6 — >■ cud 
mode, V cs and the c and s quark fields are replaced by V u d and the w and 
<i quark fileds, respectively. It is simpler to work with the operators 0± = 
|(0 2 ±Oi) and their corresponding coefficients c± = c 2 (/x) ±Ci(/i), since they 
are multiplicatively renormalized. The expressions of the Wilson coefficients 
c± are referred to 

To simplify the analysis, we route the transverse momentum p±_ of the b 
quark through the outgoing c quark as in the semileptonic cases, such that 
the c quark momentum is the same as before. For kinematics, we make 
the correspondence with the c (u) quark carrying the momentum of the 
massive (light) lepton r (e and fi) and with the s and d quarks carrying the 
momentum of v. The scaling variables are then defined exactly by Eq. (Q). 
The factorization formulas for the nonleptonic decay rates are quoted as || 

r NL mi /•,,,/•! , „ riv c +i . n c -i 



r 1 POO 

dxdydy / dz bdb 



r 2tt 

xf B (z)J c (x, y, y 0l z, b)H(x, y, y , z) 



x exp 



(10) 



The B meson distribution function being universal, is the same as that 
employed for the semileptonic decays. 

The double logarithms contained in a light quark jet have been resummed 
into the Sudakov factor exp[— 2s(pj, b)} 0, j — s for the b — > ccs mode and 
j — u, d for the b —> cud mode, where the upper bound pj = (pj + pj) of the 
Sudakov evolution is chosen as the sum of the longitudinal components of pj. 
We refer to [IT] for the expression of the exponent s. In the numerical analysis 
below, exp (—2s) is set to unity as pj < 1/6. In this region the outgoing quarks 
are regarded as being highly off-shell, and no double logarithms are associated 
with them. These quark lines should be absorbed into the hard subamplitude, 
instead of into the jet functions. The Sudakov factor, though important only 
in the end-point region which is not essential for the calculation of total decay 
rates, improves the applicability of PQCD Kued below. 

Another exponential factor comes from the RG summation of the single 
logarithms \n(tb), which is written as ^ 



Ji/b fx 



(11) 



7j = —2a s /7T and 75 = —Cpa s /Tr being the anomalous dimensions of the 
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light-quark jets and of the distribution function, respectively. The summa- 
tion of the single logarithms \n(M w /t) leads to the Wilson evolution denoted 
by c±(t). The scale t is chosen as the maximal relevant scales, 

r 



t = max 



Pv 



(12) 



Note that t, depending on the quark kinematics, probes the specific dynamics 
of heavy hadron decays. That is, PQCD effects differ between the b — > cud 
and b — > ccs decays and between the B meson and A& baryon decays. The 
above single-logarithm evolutions play an important role in enhancing the 
nonleptonic branching ratios and thus in lowering the semileptonic branching 
ratios. This observation is basically consistent with that in Since the 
large b region is Sudakov suppressed by exp(— 2s) as stated before, t remains 
as a hard scale, and the perturbative evaluation of the hard subamplitude H 
is reliable. 

Without the PQCD evolutions, our formalism is equivalent to the HQET- 
based OPE. This is obvious by expanding the distribution function in 
terms of its moments in the factorization formulas for the semileptonic 
decays. though constructed on nonlocal oparators, is constrained up 
to the second moment, the same as in the HQET approach. The large 
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logarithmic radiative corrections have been summed to all orders into the 
evolution factors as shown in Eq. ([II]). The logarithms ln(M&/M c ) in the 
semileptonic decays are less important and neglected as in Eq. @. The hard 
subamplitudes H are the initial conditions of the single-logarithm evolutions. 
In the HQET approach radiative corrections are evaluated exactly to finite 
orders. Hence, the difference between the two approaches resides only in 
the treatment of PQCD effects. We shall show that PQCD effects, being 
pro cess- dependent, are responsible for the explanation of the lifetime ratio 
r(A b )/r(B d ). 

The above formalism can be generalized to the A& baryon decays straight- 
forwardly, for which the hard subamplitudes H, the final-state cut J c , and the 
Wilson coefficients c± remain the same but with the kinematic variables of 
the A;, baryon inserted. For the A& baryon distribution function, we assume 
the similar parametrization 

^w= V-^+V (13) 

The moment associated with f\ b (z) need to be treated as a free parameter, 
since spectra of A& baryon decays are still absent, and f\ b has not been 
determined yet. The relations of the parameters N, a and e to the first three 
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moments of f\ b {z) are demanded by the HQET-based OPE: 



i h A(A 



q f Ab (z)dz = 1 , J o (1 - = ^ + °( A QCdM 2 



\l-z? 



2 h b (z)dz = ^ (A 2 (A b ) - ±A X (A*)) + (9(A 3 QCD /Mij, (14) 

with Ai(Aft) = (A 6 |6(?D ± ) 2 6|A 6 ) and 

M A6 = M b + A(A 6 ) - . (15) 

Note that Ai and A are not physically well defined quantities due to the 
presence of infrared renormalon ambiguity. They may bear very different 
values in different schemes. 

We argue that only Ai(A&) is a free parameter. The moments of the B 
meson distribution function obey a mass relation similar to Eq. (|i~5|), 



M B = M b + A(B)-^. (16) 

Since the parameters A(B) = 0.65 GeV and X\(B) = —0.71 GeV 2 have been 
extracted from the available B — > Xgj spectrum (see Eq. (§)) and Mb = 
5.279 GeV is known, M b = 4.551 GeV is fixed by Eq. flT|). Substituting 



this M b into Eq. (fL|) with M Ab = 5.621 GeV, A(A 6 ) will be derived, if 
Ai(Afe) is chosen. Using the values of A(A&) and Ai(A&), combined with the 
normalization of f\ b (z), we obtain the corresponding parameters N, a and e. 
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Several reasonable choices of Ai(Aft) are made, and results for M c = 1.6 
GeV, M D = 1.869 GeV and M T = 1.7771 GeV are obtained. Our predictions 
are not sensitive to the change of M c . The difference of the results for M c = 
1.5 GeV from those for M c = 1.6 GeV is less than 5%. We adopt the leading- 
order Wilson coefficients with the number of active quarks rif = 5. On the 
experimental side, the CLEO group reports Bsh = (10.19 ± 0.37)% and 
(n c ) = (1.12 ± 0.05) p[, while the CERN e+e" collider LEP measurements 
give B SL = (11.12 ± 0.20)% and (n c ) = (1.20 ± 0.07) jLj. The lifetimes of 
the Bd meson and of the A& baryon are r(Bd) = (1.56 ± 0.06) x 10~ 12 s and 
r(Ab) = (1.23 ±0.05) x 10 -12 s, respectively, from the LEP measurements [[]]. 
Recent CDF result yields (1.32 ± 0.16) x 10 -12 s for the A b baryon lifetime 
PH . The lifetime ratios T(A b )/r(B d ) are then 0.79 ± 0.06 from LEP and 
0.85 ±0.11 from CDF. 

The conventional values of Ai obtained from QCD sum rules are A™ eson ~ 
A bar y on = ^ ± Q ^ GeV 2 _ Choosing Xl (A b ) = -0.4 GeV 2 (with the 

correspondiong A(A^) = 1.03 GeV), we derive the lifetime ratio about 0.81, 
which is much smaller than that from the OPE approach. If setting Ai(A&) 
to Ai(.B) = -0.71 GeV 2 (with the correspondiong A(A fe ) = 0.99 GeV) that 
was extracted from the B —>■ X s ^ decay, we arrive at a ratio 0.78 almost 
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the same as the central value of the LEP result. If further increasing the 
magnitude of Ai(A b ) to —1.1 GeV 2 (with the correspondiong A(A fe ) = 0.95 
GeV), an even lower ratio about 0.75 is reached. From the above analyses in 
a wide range of the HQET parameters, it is obvious that the lower lifetime 
ratio is mainly due to the PQCD effects. Our predictions for the absolute 
lifetimes of the B meson and of the A b baryon prefer a larger \V c b\. Adopting 
\V cb \ = 0.044, we derive r(B d ) = 1.56 x 10~ 12 s and r(A 6 ) = 1.26 x 10~ 12 
s, 1.22 x 10~ 12 s, and 1.18 x 10~ 12 s for Ai(A b ) = -0.4 GeV 2 , -0.71 GeV 2 , 
and —1.1 GeV 2 , respectively, which are consistent with the data. We have 
adopted |V CS | = \V u d\ = 1.0 in the present work. 

In Table I we compare the branching ratio of each mode in the B meson 
and Ab baryon decays in terms of the quantities r TV = BR{b — > ctu) / BR{b — > 
cZP), r ud = BR(b -> cud)/BR(b -> dv), and r cs = BR(b -> ccs)/BR(b -> 
cW), and observe that B^l, (n c ) and all r's except r u d are enhanced a bit. It 
is interesting to examine these tendencies in future experiments. The results 
listed in Table I are insensitive to the variation of the distribution functions 
(less than 10% under 40% variation of the moments), confirming that the 
branching ratios are mainly determined by the PQCD effects. We emphasize 
that the predictions for all the branching ratios of the B meson decays are 
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consistent with the experiment data |16| . 
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Table I. Comparision of t tv , r u d, r cs , B^l, and (n c ) in the B meson and 
A fe baryon decays for various Ai(A 6 ). The lifetime ratios r(A b )/r(_B rf ) are also 



listed. 







7~tv 


Tud 


7~cs 




K> 


T(A b )/r(B d ) 


MB) = 


-0.71 GeV 2 


0.224 


5.98 


1.64 


10.16 


1.17 




Ai(A 6 ) = 


-0.40 GeV 2 


0.254 


5.14 


1.75 


10.94 


1.19 


0.81 


Ai(Aft) = 


-0.71 GeV 2 


0.261 


5.26 


1.68 


10.87 


1.18 


0.78 


Ai(A 6 ) = 


-1.10 GeV 2 


0.266 


5.36 


1.68 


10.76 


1.18 


0.75 
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